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Abstract 



The post Gaussian effective potential in D = 3 dimensions and the Gaus- 
sian effective potential in D = 2 + e are evaluated for the Ginzburg-Landau 
theory of superconductivity. It is shown that, the next order correction to 
the Gaussian approximation of the Ginzburg-Landau parameter n is signifi- 
cant, whereas contribution from the two dimensionality is rather small. This 
strongly indicates that strong correlation plays a more dominant role than 
the two dimensionality does in high T c superconductivity. 
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I. INTRODUCTION 



The Ginzburg-Landau (GL) theory of superconductivity [1] was proposed long before 
the famous BCS microscopic theory of superconductivity was discovered. A few years after 
the appearance of the BCS theory, Gorkov derived [2] the GL theory from the BCS theory. 
Amazingly , the GL theory has played a significant role in understanding superconductivity 
up to now. It is highly relevant for the description of high- T c superconductors, even though 
the original BCS theory is inadequate to treat these materials. The success of the GL theory 
in the study of modern problems of superconductivity lies on its universal effective character 
in which the details of the microscopic model are unimportant. 

Even in the level of meanfield approximation (MFA), the GL theory gives significant 
information such as penetration depth (I) and coherence length (£) of the superconducting 
samples. Many unconventional properties of superconductivity connected with the break 
down of the simple MFA has been studied both analytically [3] and numerically using the 
GL theory [4]. Particularly, the fluctuations of the gauge field were studied recently by 
Camarda et. al. [5] and Abreu et. al. [6] in the Gaussian approximation of the field theory. 
The effective mass parameters of the Gaussian effective potential (GEP), Q and A , were 
interpreted as inverses of the coherence length £ = 1/Q and of the penetration depth £ = 
1/A, respectively. 

In this note, we take one step further estimating corrections to the Gaussian effective 
potential for the U(l) scalar electrodynamics where it represents the standard static GL 
effective model of superconductivity. Although it was found that, in the covariant pure 
\cf) A theory in 3 + 1 dimensions, corrections to the GEP are not large [7] , we do not expect 
them to be negligible in three dimensions for high T c superconductivity, where the system 
is strongly correlated. 

Apart from the strong correlation , another important factor, which one should consider 
for high T c superconductivity, is the dimensionality of the system. It is well known that, most 
of the high T c superconducting materials have layered structures, which strongly suggests 
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two-dimensional nature of high T c superconductivity. In order to test relative importance 
of the dimensionality contribution compared to the post Gaussian corrections, we shall also 
study the case of fractal dimension, D = 2 + e. 

The paper is organized as follows: in Section II the GL action is introduced and the 
post Gaussian approximation is applied; in Section III, the theoretical results for D = 3 
and D = 2 + e will be compared to existing high T c experimental data. The results are 
summarised in Section IV. 

II. POST GAUSSIAN EFFECTIVE POTENTIAL IN D=3 DIMENSION 

We start with the Hamiltonian of the GL model in Euclidean D-dimensional space given 



where ip and A are the complex scalar and the static electromagnetic fields, respectively; 
m, A and e are the input parameters of the model. 1 We introduce natural units employing 
£o (coherence length at zero temperature) and T c as length and energy scale, respectively, 
through the transformations : 



by [8] 



+ imV + A^ (3 - D) H 4 } , 



(2.1) 



m — y m£, 



— l 
o > 



■o X ) 



X -> X1JQ, 



(2.2) 



Eq. (2.1) is now rewritten as, 




(2.3) 



i 



jjl is introduced to make A and e dimensionless. 
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In accordance with refs. [6,5], we apply tranverse unitary gauge and express the partition 
function as 

Z = J V(f)VA T exp{- J d D xH + J d D xj<f) + (j A A)} (2.4) 
where the Hamiltonian density is 2 

H = i(V x Af + 1(V0) 2 + ^m 2 2 + A0 4 + ^e 2 2 A 2 

+ 1(V/) 2 . (2.5) 

We have introduced a gauge fixing term, with the limit e — > being taken after the calcula- 
tions are carried out. In Eq. (2.5) A stands for the transverse gauge field and ip = (fiexp(ij) . 

To obtain the free energy density, V q q = T fV (effective potential), we introduce a shifted 
field — > + 0o an d split the Hamiltonian into two parts: 

H = H + H int , (2.6) 

— # 

where H is the sum of two free field terms describing a vector field A with mass A and a 
real scalar field with mass flo: 

iJ = i(Vxl) 2 + iA 2 l 2 + l(Vl) 2 

+ i(V0) 2 + ^ 2 2 . (2.7) 
The interaction term then reads 

F int (0, A) = ]T - 1a 2 ! 2 + ie 2 i 2 (0 + O ) 2 , (2.8) 

n=0 

where 

v = im 2 2 + A^ 3 "^ , Vl = m 2 o + 4\^~ D )$ , 

U2 = I( m 2-ng) + 6A M ( 3 - D V§ ) ^3 = 4A^ 3 - D )0 O , (2-9) 

V4 = . 



2 From now on, we denote A/t/ 3 D ) and e 2 ^ 3 D ^ as A and e 2 , respectively, for simplicity. 
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Now performing explicit Gaussian integration in Eq. (2.4), one obtains 

Z = exp{- J d D xH- mt (5/Sj,S/Sj A )} J V4>VA^ V { 
- J d D xH +j ( f } + j A A)} 

= [det J D - 1 ]-5[detG 1 ]-^ex P {- J d D xH [nt (5/6j,5/5j A )} 

x exp{jD j/2}exp{j A G j A /2} , (2.10) 

where in momentum space 

D (p) = l/(p 2 + nt), G (p) = 2/(p 2 + Al). (2.11) 

To calculate partition function in post Gaussian approximation, we use the method intro- 
duced in refs. [9,10] and introduce the so called primed derivatives: 



*5j{x)> x 5j(x)> K 5j A {x)> x 8j A (x)' 



5 2 - 5 2 

( f . 9 , J = = f .o/ x - D 0(X ,X) 



( ^R )/ = ^ 2) = ^R- Go(x ' x) ' 



( ^y.^ ) = _|__3 Go(x ^ ) 



(5j 4 (x) dj 4 (a;) 5j J (x) 



t)^ 

: w ) ' Efl = «w" 6C,(l,l) «w 



where = ^ d D yD (x ,y)j(y) and itU(z) = ^ d D yG (x ,y)j A (y), so that 



ii n) exp{j J D j72} = R n (x)eMjD j/2} , 

(2.13) 

Bi n) exp{j A G j A /2} = R n A (x)exp{j A G j A /2} . 
Now it can be shown that [9,10], the Gaussian part of Z can easily be isolated as follows: 

Z = Z G AZ , 

Z G = exp{-/ 1 (fi) - i/i(A) -vo- v 2 I (Q) + 3v 4 I 2 (Q) 

+ (A2 + e 2 / (fi)-e^)/ (A)}, 

AZ = exp{-t; 2 i( 2 ) - v 3 A^ - v 4 i (4) 
1 



2 



e 



2 ^ 



A 2 )£^ - e 2 o # (2) ^ (1) (2.14) 



ie 2 J B( 2 )i( 2 )}exp{j J D J 72}exp {j a Gj a /2} , 



where 



£>(p) = l/(p 2 + n 2 ), G(p) = 2/(p 2 + A 2 ), 

= n§ + 12v 4 I (n) + 2e 2 / (A), A 2 = Ag + e 2 / (ft) 
In the above, following integrals are introduced 

Io(M) 1 



(2.15) 



(2ir) D (M 2 + p 2 ) ' 



(2.16) 



(2^ 

From Eqs. (2.14), (2.15) and (2.9), one gets the following Gaussian effective potential: 

V G = - In Z G = h(n) + i/x(A) + v + v 2 I {Q) 
- 3vJ*(n) - (A 2 + e 2 / (fi) - e 2 2 )/ o (A) 

= h(Q) + h ± (A) + im 2 2 + \4 (2.17) 

+ i/ (fi)[m 2 - Q 2 + 6A/ (fi) + 12 A0 2 ,] 

+ Jo(A)[-A 2 + e 2 / o (fi) + e 2 2 ]. 

Note that, the last equation is exactly the same as it is in refs. [5,6]. The post Gaussian 
effective potential 
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V eG = V G + AV G (2.18) 

includes a correction part AVq : 

AVb = -InAZ = -\n{exp[-5W] 

x exp{jDj/2}exp {j'aGjV 2 }I:;=o,m=o} 

= _ ln{l - 5Wexp{jDj/2}ex P {j a Gj a /2}\ 1=0ja=0 
5 2 W 2 

+ -y-exp{jDj/2}exp {j a G Ja /2}\ j=0>Ja=0 + ...)} 

= 5AV^ (B) + 5 2 AV^ 2) (B) + . . . , (2. 19) 

^ = V2 A^ + ^(3) + W4 iW + I( e 2 2 _ A 2 )j g(2) 
+ e 2 0Oj B(2)l + 1^5(2)^(2) 

dj 2 

Here we have introduced an auxiliary expansion parameter 5 to be set equal to unity after 
calculations. 

The first order term AVq\b) in this equation will not contribute to the effective poten- 
tial, i.e., AVq\b) = 0, due to the relations (2.13). The next term of order 5 2 gives the first 
nontrivial contribution to the post Gaussian effective potential. The explicit calculations 
give 

AV G = [~ e 4 1 2 (A) - 18 1 2 (Q) A 2 ] O 4 + {-3 A I 2 (Q) 
x [-n 2 + m 2 + 2 Iq(A) e 2 + 12 A I (^)] 

- e 2 I 2 (A) [-A 2 + e 2 I (Q)] - 8 A 2 l 3 (Q, Q) (2.20) 

- \ e 4 1 3 (A, fi)}0 o 2 - \ I 2 (fi) [-{I 2 + m 2 + 2 1 (A) e 2 

o o 

+ 12 A I (fi)] 2 - i I 2 (A) [-A 2 + e 2 1 (fi)] 2 

-l e 4 i 4 (A, n)-^x 2 un,n), 

where following loop integrals were introduced, 
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(2.21) 



I 4 (M 1 ,M 2 ) = 



1 f d D kd D pd D q 



(2tt) 3D J (k 2 + M?)(p 2 + Ml)(q 2 + Mi) 



1 



x 



((k + p + q) 2 + M 2 ) ' 



For D = 3 — 2e, these integrals were calculated in dimensional regulariztion in ref. [11]. For 
completeness explicit expressions are given in the Appendix. The appropriate counter terms 
to cancel the divergences coming from the integerals are also presented in the Appendix. 
The parameters Q and A are determined by the principle of minimal sensitivity (PMS) : 



= 576 A vr 2 A 2 $(2 A + ti ) + 8ir$j{-12 A e 2 A 3 

+ [-72 7rA^ 2 + (-6Ae 2 + 192A 2 )^ + 247rm 2 A]A 2 

+ A [-36 7rAf] 3 + (10e 4 + 96A 2 )^ 2 

+ 12 tt^ Am 2 ] +n 3 e 4 } +20 A H 5 vr 2 

+ (48^ A A +40A 2 tt 2 -2e 4 )tt 4 

- 4Af2 3 [-24vr A A - A vr e 2 + 2 e 4 ln2 + e 4 




+ 67r 2 m 2 -24A 2 ln^- 



(2.22) 



- 8 A 2 n 2 [- A vr e 2 + e 4 ln^ + 6 n 2 m 2 

- 2 e 4 ln /X ^ 2 - NO - 15 A 2 + 2 e 4 ln2 



+ 



+ 96 A 2 In2-24A 2 ln-^1 

+ A(2A +fi)(A e 2 -2 7rm 2 ) 2 = 0; 
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= 16 fi n 2 e 4 ${n + 2 A) (2.23) 

+ 8e 2 7r0§{(-8fi 7r + 12 A) A 3 

+ 2U A 2 (-2n vr + 4e 2 + 3A) - e 2 2 (fi + 6 A )} 

+ e 4 4 + 2An 3 e 2 (2^A+5e 2 ) 

+ 4A 2 ^ 2 [4vr 2 A 2 + 2vre 2 A+6e 4 

+ 2e4 ' n (AW- 4e4 " ,2 - 3e2A1 

+ 2A 2 [16^ 2 A 3 + 7e 4 A 

+ 8A e 4 rn - ^ - x9 - 12 A e 2 A-16A e 4 ln2 
(A 

+ 2e 2 m 2 Tr]n +4A 3 e 2 (-A e 2 + 2 7rm 2 ) = 0, 



where we denote optimal values of Q and A by Q and A , respectively, and 0o is a stationary 
point defined from the equation: 



{(-47re 4 + 647r 2 A A) Q - 144vrA 2 A}0, 



o 



d(f> 

+ (e 4 - 367rAA)fi 2 + 2fi A[-24A 2 ln-^ 

a 2 

+ 48 A 2 ln3 - 2 A 7T e 2 + e 4 ln^ - 2 e 4 ln2m 2 

A 1 

+ 8l2 - 2e4 '"(2AT!TF- 6A21 (224) 
+ 6AA (A e 2 - 2vrm 2 ) = 0. 

Note that, in the Gaussian approximation, the gap equations (2.24) are reduced to simple 
forms: 

BV 

- 7 ^ = fi 2 7r-6fiA-Ae 2 + 2m 2 7r = 0, 
oil 

(2.25) 

^ = 8 A 2 7T A - e 4 A - Ail e 2 A + 2e 2 m 2 tt = . 
oil 



III. COMPARISION WITH EXPERIMENTAL DATA FOR D = 3 AND D = 2 + e . 

The solutions of the Eqs. (2.24) are related to the experimentally measured GL parameter 
k as k — £/£ — Cl /A . We make an attempt to reproduce recent experimental data on k(T) 
[12] for high -T c cuprate superconductor Tl 2 Ca 2 Ba 2 Cu 3 O w (T£ - 2223). 

For this purpose, we adopt usual linear T dependence of parametrization of m and A as: 

m 2 = mg(l — t) + rm 2 , 

A = A (l-r)+rA c , (3.1) 
r = T/T c , 

and calculate k by solving nonlinear equations (2.24) or (2.26). Due to the parametrization 
(3.1), the model has in general five input parameters: m 2 ,, A , ra 2 , A c and e. The last 
parameter is related directly to the electron charge: e 2 = 16nakBT c !; /hc, where a = 1/137, 
£o is a coherence length at T = , and T c the critical temperature. The experimental values 
for the cuprate T£ — 2223 are £ = 1.36nm and T c = 121.5K. The parameters wig and A are 
fitted to the expeimental values of £ and £ at zero temperature: £o = 1.36nm, £ = 163nm. 
In dimensionless units, (2.2), we have Q = Q(t = 0) = 1 and A = A (r = 0) = £oAo = 
0.0083 which are used to calculate m 2 , and Ao from coupled equations (2.24) (or (2.26) in 
the Gaussian case). The parameters m 2 . and A c are fixed in the similar way. Actually the 
quantum fluctuations shift m 2 from its zero value given by MFA. On the other hand, the 
exact experimental values of m 2 and A c are unknown, since the GL parameter at T = T c is 
poorly determined. For this reason, we used the experimental values of £ c and £ c at very close 
points to the critical temperature taking r c = 0.98 which gives Cl c = f2 (r c ) = l/£ c = 0.128 
and A c = A(t c ) = 1/4 = 0.0043. Then solving the equations (2.24) (or (2.26) in the 
Gaussian case) with respect to m c and A c , we fix the input parameters. Their values for the 
Gaussian and the post Gaussian cases for D=3 are summarized in Table 1 3 . 



3 All parameters are given in dimensionless units. See Eq. (2.2). 
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After having fixed the input parameters, the temperature dependence of Cl(r), A (r) as 
well as the GL parameter k = Cl (r)/A (t) are established by solving the gap equations (2.25) 
and (2.24) numerically for the Gaussian and the post Gaussian approximations, respectively. 
The results are presented in Fig.l , where solid curve corresponds to the post Gaussian and 
dotted one to the Gaussian approximation. It is seen from the figure that corrections to the 
Gaussian approximation are significant, and in the right direction, although the discrepancy 
from the experimental values is still substantial. 

On the other hand, a better agreement with the experiment has been obtained even on 
the level of the Gaussian approximation by the authors of ref. [5]. However, they introduced 
a cut off parameter A as a characteristic energy scale of the sample to make the divergent 
integrals I and I\ finite. We beleive that, the better agreement is a result of introducing 
this rather arbitrary additional parameter. It should be noted that, in the present approach, 
there is no such additional adjustable parameter. Here we used dimensional regularization 
in which we put /i — Cl . It was found that the behavior of k{t) does not depend on /i: 
Another value of /i, e.g. \i = 2Q leads to another set of input parameters {ml,m c , A , A }, 
but to the same behavior for k(t). 

Clearly, the solutions of nonlinear gap equations are not unique. In numerical calculations 
we separated the physical solutions by observing the sign of 4>q and that of the effective 
potential at the stationary point : V e g(0 o )- The temperature dependence of these two 
quantities are presented in Fig. 2. It is seen that 4>1 (solid line) is positive in the large range 
of r and goes to zero when r is close to r = 1. Similarly, the depth of the effective potential 



TABLE I. Input parameters of the GL model. 





ml 


Ao 


o 


A c 


Gaussian 


-0.456 


0.046 


0.0013 


0.002 


Post. Gaussian 


-0.525 


0.050 


0.0017 


0.008 
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at the stationary point, V^g(0o), becomes shallow when r — > 1 and vanishes at T = T c . 

All the above numerical calculations were made in D=3 dimension. On the other hand it 
is widely known that, most of high T c cuprates have layered structures with 2D CuOi planes 
which play an essential role in the high T c superconductivity. Therefore, it is nessesary to 
consider the dimensional contribution in the calculation so that relative importance between 
the post Gaussian corrections and the two -dimensional character can be assessed. For this 
purpose, we consider the case of D = 2 + 2e (e ^ 0) in the Gaussian approximation. The 
effective potential is given by the Eqs. (2.16) and (2.17) where the integrals are explicitely 
written as: 



U M ) = -— {[— vr 2 + - In 2 e - In ^- + -} , 
0V ; 4tt 1[ 12 2 M 2J M 2 e 



(3.2) 



Taking derivatives of Vq and using (3.2) leads to the following gap equations: 
dV G] 



on 



\(dV G /d<f> )=0 



6\ + e z -£[ttQ 2 + 2nm z 
+ (e 2 + 12A)ln^+e 2 ln|l] 

+ O(e 2 )=0, (3.3) 



' ° \(dv G /d<f> )=o = 4e 2 A + e 4 - 2e[AA 2 ttA + e 2 m 2 Tc 
+ e 2 (e 2 + 2A)ln|^ 

+ 2Ae 2 ln^]+O(£ 2 ) = 0. (3.4) 

The parameters ro§, A , m 2 and A c in Eqs. (3.1) - (3.4) were adjusted to their experi- 
mental values in the same way as in the previous case. As a result we obtain: 



dA 
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1.106 - 22.827 e 
1.106 + 36.443 e' 



1.106 + 36.443e 



0.869 e 




-511.839 e 4 + 15.368 g 3 + 0.708 e 2 - 0.00645 e + 0.000016 
140.569 e 2 + 2982.92 e 3 + 1.9283 e + 14677.46 e 4 



(3.5) 



0.202 e 3 + 0.0042 e 2 - 0.00005 e 



1.405 e + 29.829 e 2 + 0.0192 + 146.77 e 3 ' 



In Fig. 3 and Fig.4, we present m 2 and A vs. e, respectively, given by Eqs. (3.1) and 
(3.5). One notes that for small values of e (0 < e < 0.048) m 2 becomes positive . Bearing 
in mind that, in the GL model the phase transition occurs where m 2 changes sign (or more 
exactly the superconductive phase holds only for m 2 < 0), it shows that, in the present 
approximation scheme, there is no phase transition in D = 2 + 2e dimension for very small 
values of e. This smallness of the e value indicates reliability of the present post Gaussian 
approximation. Note that , A remains positive on the whole range of e (Fig.4). 

The GL parameter «(r) = 0(r)/A(r) given by Eqs. (3.1) - (3.5) in D = 2 + 2e 
case (e = 0.1) is plotted in Fig.l (dotted line). We find that, surprisingly, the two curves 
from the Gaussian approximation almost coincide, thus indicating that the dimensionality 
contribution is not significant as far as the GL parameter is concerned. 



In the present article we have carried out calculations on the Ginzburg-Landau effective 
potential beyond the Gaussian approximation. The result is used to obtain the Ginzburg- 
Landau parameter, k, and compared with existing high T c superconductivity data. It was 
shown that the post Gaussian correction which is believed to originate from strong corre- 
lation is substantial. In order to estimate the contribution from the two dimensionality of 
high T c superconducting materials, we have carried out calculations for D = 2 + 2e in the 
Gaussian approximation. The result shows that the dimensionality correction to the three 
dimensional Gaussian result is rather small, although there remains possibility that a post 



IV. SUMMARY 



13 



Gaussian correction at D = 2 + 2e is much larger than that at D = 3, thus making the 
theory closer to experiment. This remains as a future study. 
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Appendix 

A. Explicit expression for divergent integrals. 

Here, we bring explicit expressions for the divergent integrals defined as: 

/o(M) = -^[l + £ (2 + ln ^)+0(. 2 )], 



M 3 u 2 



(A.l) 



+ ln-^ -6m2 + - + 0( £ )] 
Mi e 

+ Mi ^ M 2 } . 

In D=3 dimension the integrals J 3 and J 4 are divergent. Following counter terms were 
introduced to the Hamiltonian [9,13,14]: 

H cont = B</) 2 m 2 /2 + CA0 4 + De 2 </) 2 A 2 (A.2) 

+ EA 2 A 2 (l-5)/2, 



where in the minimal subtraction (MS) scheme 

B = B 1 5 + B 2 5 2 + . . . , C = d6 + C 2 5 2 + . . . , 
D = DiS + £> 2 <5 2 + . . . , E = E 1 5 + E 2 5 2 + ... , 

B L = 0, Ci = 0, L>i = 0, = , 

3A 2 e 2 e 4 
-B2 = : o , Co = 0, Z?2 = — : , E\ 



(A.3) 



27r 2 e ' " ~ Aire ' 167r 2 e 
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FIG. 1. The GL parameter, k, in the Gaussian approximation in D = 3 (the dotted line) and 
D = 2 + e (the dashed line) cases. The solid line represents the PostGaussian approximation for 
D = 3 case. 
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FIG. 3. The parameter m 2 of th GL model given by Eqs. (3.1) and (3.5) for D = 2 + 2s. 
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FIG. 4. The same as in Fig. 3 but for A 
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